Generalized Haldane Equation and Fluctuation Theorem in the Steady State Cycle 

Kinetics of Single Enzymes 
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Enyzme kinetics are cyclic. We study a Markov renewal process model of single-enzyme turnover 
in nonequilibrium steady-state (NESS) with sustained concentrations for substrates and products. 
We show that the forward and backward cycle times have idential non-exponential distributions: 
Q+(t) = @-(t). This equation generalizes the Haldane relation in reversible enzyme kinetics. In 
terms of the probabilities for the forward (p+) and backward (p_) cycles, fesTln(p+/p_) is shown to 
be the chemical driving force of the NESS, A/Lt. More interestingly, the moment generating function 
of the stochastic number of substrate cycle v(t), (e~ Xv ^) follows the fluctuation theorem in the 
form of Kurchan-Lebowitz-Spohn-type symmetry. When A = A/i/fcsT, we obtain the Jarzynski- 
Hatano-Sasa-type equality: (e~' / ' t ' AA1 /' :sT ) = 1 for all t, where i/A/i is the fluctuating chemical work 
done for sustaining the NESS. This theory suggests possible methods to experimentally determine 
the nonequilibrium driving force in situ from turnover data via single-molecule enzymology. 
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Most biochemical reactions in a living cell have nonzero 
flux J and nonzero chemical driving force A[i. The 
nonequilibrium state of such a reaction is sustained by 
continuous material and energy exchange with and heat 
dissipation into its environment 0. To understand the 
state of a biochemical network in an open environment, 
hence, it is necessary to be able to experimentally mea- 
sure both J and A/i in situ. A large literature exist on 
measuring J, but none exists on directly measuring Afx. 
One could in principle compute A/i from in situ measure- 
ments of the concentrations of the substrate and product 
of a reaction if its equilibrium constant is known Q . Al- 
ternatively, one should be able to obtain A/x from fluctu- 
ating cycle kinetics of a single enzyme directly. This pos- 
sibility has been recently investigated in term of stochas- 
tic simulations Here we exam this idea through a 
novel analytical model. 

Enzyme kinetics are complex mainly due to the many 
possible intermediates in the form of enzyme-substrate 
complexes. Recent laboratory measurements with high 
resolution at the single-molecule levelgive the waiting 
time distributions for enzyme cycles Q. This motived 
the present Markov renewal process (MRP) model, also 
known as extended kinetics model in the theory of motor 
proteins [5| • In terms of the MRP, the kinetics of a single 
enzyme becomes a stochastic sequences of forward and 
backward cycles as function of time. We shall denote 
the number of forward and backward cycles by v+ (t) and 
V-(t), as shown in Fig. 1. 
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FIG. 1: The solid line illustrates ideal data on single enzyme 
cycling as function of time, v(t), which can be decomposed 
into v+(i) and V-(i), shown in dashed and dotted lines. The 
starting positions are arbitary. AT+ and AT- are forward 
and backward cycle times. 



It is obvious that the cycle time distributions give in- 
formation on the kinetics. In this letter we show that the 
key nonequilibrium thermodynamic quantity, A/i, can 
be obtained from stochastic data on single-enzyme cy- 
cle v(t) = v+(t) — V—(t) via two novel equalities 



A/i = fc fl Tln((i/+(t)>/(i/_(t)» 
( e -"(*)A/i/fc fl Tj = ^ vtj 



(1) 
(2) 



2 



where (■ ■ • ) is ensemble average for repeated measure- 
ment of v(t) in steady-state. Eq. generalizes a result 
well known for one-step chemical reactions 0,0- Eq. (J2J 
is a version of the fluctuation theorem (FT) in nonequilib- 
rium statistical mechanics. The FT for the probability 
distribution of entropy production of a noncquilibrium 
steady-state (NESS) was first discovered in deterministic 
dynamical systems []|. Kurchan, Lebowitz and Spohn 
(KLS) introduced a parallel theory in terms of stochastic 
dynamics || which is more appropriate for single enzyme 
experiments 0, 0, 0] . It was shown that the generating 
function, i.e., an exponential average, of a work func- 
tional W(t) possesses a certain symmetry in the limit of 
t — > oo. Crooks introduced a heat functional Q(t) and 
showed the similar symmetry is valid for all finite t [Tlj : 



is a key result of this work. It has to do with micro- 
scopic reversibility. There are experimental edvidcncc 
for it, as well as theoretical models proving equal mean 
time (AT+) = (ATI) [H E^. We shall give a proof 
for the equal distribution with sequential enzyme kinet- 
ics. The proof for more general systems will be published 
elsewhere |l8j . 

The Markov renewal process (MRP) model. — Detailed 
kinetic scheme of an enzyme catalyzed biochemical reac- 
tion A ^ B is usually very complex . But if one con- 
cerns only with the net number of steady-state turnovers 
from A to B, v(t), it can be represented by a continuous- 
time, discrete-state one dimensional random walk with 
cumulative cycle time distribution functions Q±(t) for 
the forward and the backward stochastic transition times 



cx(t) = a-x(t) where c x {t) = ( e -AQ(*)/*aiy Since at+ and AT _. Q± ^ = ^ e± ( 00 j _ ^ and 0± ( t ) 



and W(t) differ by a stationary term while both increase 
without bound, Crooks' result immediately yields that of 
KLS's. The symmetry in the generating function implies 
FT for Q 0. 

The symmetry implies that ln(e -< 20/' CBT ) = 0. This 
is analogous to the Jarzynski equality which is sur- 
prising since (Q(t)) = — lne~^^'' kBT is the mean heat 
dissipated from the NESS which certainly ^ 0; it should 
always > 0. The Jarzynski equality provides the possibil- 
ity to obtain a function of state such as free energy from 
an nonstationary heat functional Q(t) with finite t. This 
was proposed and experimentally tested for mechanical 
work functional on single biological macromolecules such 
as RNA [HQ- 

The difference between the FTs for Wit) in the limit 
of infinite t and for Q(t) with any finite t is crucial to 
real experiments. In heuristic thermodynamic terms, the 
work functional Wit) is related to the Afj,° of a reac- 
tion and the heat functional Q(t) to the A/i. While 
the former is determined by the transition rate constants, 
hence experimentally accessible in short time, the latter 
depends on the stationary probability. For cyclic enzy- 
matic turnovers, however, W — Q. Hence, the FT as- 
sociated with enzyme cycle kinetics is particularly sim- 
ple, and experimentally accessible Generalizing the 
Jarzynski equality to open-system, Hatano and Sasa's 
equality for NESS ^3 a l so suggested a possibility of com- 
puting chemical driving force for single-molecule chemical 
reactions in NESS, see flilH- 

To show Eqs. (JTJ and J2J), there are two strategies. One 
is based on traditional Markov models, i.e., master equa- 
tions, for single enzyme kinetics. Then both equations 
can be show as conseqeunces of the existing FTs 0, . 
An alternative, more insightful approach however, is to 
model the kinetics in terms of a MRP with cycle kinetics. 
In this new model, we shall show a suprizing equality be- 
tween the forward and backward cycle time distributions: 
9+(t) = 9_(t). With this equality, Eq. (QJ becomes ob- 
vious, and Eq. can be shown in elementary terms, in 
the Eqs. 0-JD) below. 

The equality 6+(r) = 9_(r) turns out to be a very 
important, unknown relation in enzyme kinetics. This 



are nondecreasing. This is a class of stochastic models 
known as MRP |2(j which has wide applications in single- 
enzyme kinetics and motor protein stepping [5|, l2lj . See 
Fig. 2 in which w±(t) = p± O ± (t) , and p+ + p— = I . p+ 
is the eventual probability of the enzyme binding 
A (B) and converting it to B (A). We shall also denote 
w(t) = w+(t) + w-(t). 

We discover that a necessary condition for the Eqs. 
Q and (J2J is the cycle-time distributions for the forward 
and backward steps being equal: Q+(t) = 0_(i). We 
call this equality generalized Haldanc equation |22| . 
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FIG. 2: (a) A schematics for an enzyme reaction converting 
substrate A to product B. In a NESS, the concentrations for 
A and B, ca and cb, are controlled through feedback by an 
experimenter. The cumulative number of B taken out by the 
time t is denoted by —oo < v(t) < oo. (b) The integer 

valued u(t) is most naturally modeled by a random walk with 
forward and backward time distributions W+(t) and W-(t) Q- 



The position of the random walker in Fig. 2b, 
v{£), models the net number of enzyme turnover. Let 
= 0, Az/i, Ai>2) Ai^, ... (Az/ = ±1) be the suc- 
cessive increments of turnover number, and T = 
0,ATi,AT 2 ,...,AT £ ,... (AT > 0) be the corresponding 
increments in time. Then the probabilistic meaning of 
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w±(t) is the joint probability for continuous AT and bi- 
nary Av: 

w±{t)=Pr{Av e = ±l,ATi<t}, (£>1). (3) 

The equation Q+(t) = @-(t) leads to w±(t) = p±w(t). 
That is the random variables Avi and ATi are statisti- 
cally independent! 

0-|- (i) = Q-{t) and microscopic reversibility. — To 
show this equality for forward and backward cycles, we 
consider a sequential enzyme reaction as shown in Fig. 
3a and a corresponding exit problem [23j | shown in 3b. 
Starting at the central position E, w + (t) and w~(t) are 
the cumulative probabilities of reaching B+E and A + E. 
Since only the first and last steps are irreversible, w+(t) 
and W-(t) both have 2n + 1 exponential terms with 
same eigenvalues, one of which is 0. Both thus can be 
written as a + a ± e^ Xlt + a 2 e~ A2t ... + a 2n e~ X2nt . With 
some straightforward algebra, it can be shown that for 
all < to < 2n HI, 



1 



d m W-(0) 



d m w+(0) 



wiw 2 ...w n 



dt" 



UlU 2 ..M n 



dt 1 ' 



(4) 



Since functions w + (t) and W-(t) are completely deter- 
mined by these initial conditions which satisfy the linear 
algebraic system, we have 



«>-(*) 



n 



(5) 



independent of t. That is 0-(i) = Q+(t). 

The meaning of the equality now becomes clear: We re- 
call that Mi and w n are pseudo-first order rate constants: 
mi = u°ca and w n = w^cb- In a chemical equilibrium: 



cb _ u\u 2 ...u n ^iu n 

CA WlW 2 ...W n -iW° 



(6) 



that is w+(t) = w-(t). Therefore in a chemical equilib- 
rium, not only on the average iu+(oc) = u>_(oo), i.e., the 
forward flux equals the backward flux, but the detailed 
kinetics for the transition time distributions has to be 
equivalent: There is absolutely no statistical difference 
between the forward and backward reactions. In a NESS 
when (0 is not held true, w+(t) ^ w-(t). But the dif- 
ference is only in the total probability p+ — w+ (oo) and 
p_ = w_(oo), the distribution functions Q+(t) — @-{t) 
still hold true. This equality is essencial to the KLS sym- 
metry below. It is known that microscopic reversibility 
has to be satisfied even when a mesoscopic system is in 
a nonequilibrium steady-state 

Kurchan-Lebowitz-Spohn (KLS) symmetry and FT.— 
For k successive number of renewal events (forward plus 
backward turnovers) within time [0, i], let us denote 
(v k ,T k ) — ~Y^l =l {Avi, ATf). The moment generating 
function for v{t) is 



g x (t) ee (e-^W) 



(7) 



= J2 e~ Xn J2Mvk = n,T k <t,T k+1 > t} 
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FIG. 3: (a) A schematic for an enzyme reaction converting A 
to B. The transition time distribution of a single enzyme con- 
verting an A to B, w+(t), and converting a B to A, w~ (t), are 
intimately related to the exit problem shown in (b) in which 
Mi and w n are pseudo-first order rate constants which depend 
on the concentrations of A and B, respectively: mi = u°ca, 
w n — w u cb- The scheme in (b) has been used to compute 
steady-state one- way flux in T.L. Hill's theory on biochemical 
cycle kinetics [a. l2a|. 



Pr{v k = n} Pr{T fe < t, T k+1 > $) 




k=0 



(p +e - x +p-e x ) Pr{T k < t.Tfc+1 > t} 



(9) 



The Eq. © is obtained due to the independence be- 
tween v k and T k . Then from Eq. © we have the KLS 
symmetry 

<?a(*)=5v-a(*), Vi, (10) 
where A* = ln(p+/p_). Furthermore, 

.g A ,(t)EE<e-^ A ^ T )= 5o (t) = l, (11) 

if ln(p + /p_) = A^/fc^T holds true. We recognize that 
v(t)A[i is the external chemical work done to the system 
in NESS. Hence Eq. i|ll|) is analogous to the Jarzynski 
equality for a cycle. 

Chemical driving force for NESS enzyme cycle. — If we 
let the t — > oo in Eq. (J5J, we have \n(p + /pJ) = A* = 
Afi/ksT, which is needed in deriving Eq. (|llfl . This 
generalizes the well-known result for single step chemical 
reactions @, |2f| to any complex enzyme reaction cycle. 

We are now also in a position to show Eq. The 
mean number of net turnovers can be computed from the 
g\(t) given in Eq. iJSJ: 



Mt)) = <!/+(*)> -<M*)> = - 



dgx(t) 



dX 



(12) 



A=0 



n= — oo fc=0 



(p + -p_)2^fcPr{T fe <t,T fc+1 >t} (13) 
— p-) x mean # of cycles in time i(14) 
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Therefore, +^4 = y-. Furthermore, in the limit of large 

t 0, {v{t)) w(p + -p_)t/(Ti>, where (T±) = tdw(t) is 
the mean time for one cycle, forward or backward. When 
p+ = p~, the steady-state flux J = lim t _» 00 (^(t))/f = 
as expected. When p + >p_, J > 0. 

Summary. — Studying enzyme-catalyzed biochemical 
reactions in situ requires methods for measuring A/i, the 
NESS chemical driving force. Currently none exists. We 
propose obtaining A/i from stochastic cycle data of a 
single enzyme molecule, v{t), via (i) an equality similar 
to that of Jarzynski and Hatano-Sasa: ( e -"W A ''/ t B T ) 
1 : or simply (ii) k B T\n{{v + {t)) / {v_ (t)}). We developed 
a MRP model for the enzyme cycles with arbitary com- 
plex mechanism, and have found an equality between the 



forward and backward cycle time distributions based on 
microscopic reversibility. This equality is a generaliza- 
tion of the what is known as Haldane relation for the 
reversible enzyme kinetics, and a recent results in 
The model enables us to establish a FT and above equal- 
ities (i) and (ii) for any t. Noting that (l/t)(v(t)) = 
J, one thus obtains both the flux J and the driving force 
A/i for a reaction in NESS from the fluctuating v(t). The 
statistical accuracies associated with these measurements 
were discussed in 
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